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The original canonical ensemble formalism for the nonextensive entropy thermostatistics is re- 
considered. It is shown that the unambiguous connection of the statistical mechanics with the 
equilibrium thermodynamics is provided if the entropic parameter l/(q — 1) is an extensive variable 
of the state. Based on a particular example of the perfect gas, it is proved that the Tsallis thermo- 
statistics meets all the requirements of equilibrium thermodynamics in the thermodynamic limit. 
In particular, the entropy of the system is extensive and the temperature is intensive. However, for 
finite systems both the Tsallis and Boltzmann-Gibbs entropies are nonextensive. The equivalence of 
the canonical and microcanonical ensembles of Tsallis thermostatistics in the thermodynamic limit 
is established. The issue associated with physical interpretation of the entropic variable is discussed 
in detail. 

PACS numbers: 24.60. Ky, 25.70. Pq; 05.70.Jk 

I. INTRODUCTION 

The main purpose of this Letter is to establish a clear way to implement the equilibrium statistical mechanics 
based on the nonadditive entropy deferent from the usual Boltzmann-Gibbs statistical one. For the first time 
this concept was formulated by Tsallis in Q . It is very well known that the conventional equilibrium statistical 
mechanics based on the Boltzmann-Gibbs entropy meets all the requirements of the equilibrium thermodynamics 
in the thermodynamic limit This is a necessary condition for self-consistent definition of any equilibrium 
statistical mechanics. In order to provide the connection of the statistical mechanics with the thermodynamics, 
the statistical entropy is usually used. From the mechanical and thermodynamical laws it allows one to determine 
a unique phase distribution function, or a statistical operator, which depends on two different sets of variables: 
the first set specifies the dynamic state of the microscopic system and the second one sets up the thermodynamic 
state of the macroscopic system. According to the Liouville and von Neumann equations, the equilibrium 
distribution function is a constant of motion which is expressed only through the first additive integrals of motion 
of the system. The ensemble averages in the statistical mechanics correspond to the concrete functions of state 
from the thermodynamics and depend only on the macroscopic variables of state. The statistical entropy as a 
function of the variables of state in the thermodynamic limit must satisfy all properties of the thermodynamic 
entropy: concavity, extensivity and so one. Note that the thermodynamic potentials of the system are functions 
fixing the norm of the phase distributions. It is known that the equilibrium thermodynamics is the theory 
defined in the thermodynamic limit. Therefore, the concept of the thermodynamic limit plays a crucial role in 
comparing the equilibrium statistical mechanics with thermodynamics. In this case, for thermodynamic systems 
the boundary effects must be neglected and only the short-range interaction forces can be taken into account. 
In the thermodynamic limit the ensemble averages with the corresponding distribution function should provide 
performance of the zero, first, second, and third laws of thermodynamics, and the principle of additivity that 
divide all variables into extensive and intensive ones. Moreover, the fundamental equation of thermodynamics, 
the Gibbs-Duhem relation, and the Euler theorem should be implemented. 

In ref. 0], the author tried to construct the equilibrium statistical mechanics based on the nonextensive 
statistical entropy. The problems arisen in the proof of the thermodynamical laws resulted in the occurrence 
of the divers variants of the Tsallis thermostatistics [sj, Q In these investigations not only the ensemble 
averages and the norm equation for the distribution function but also the Lagrange function were drastically 
redefined. Moreover, it was found that any of these variants do not satisfy the zeroth law of thermodynamics if 
the parameter q must be a universal constant. See ref. for subsequent discussions of the critique of q-entropy 
for thermal statistics. To solve this problem, in Abe et al. in the framework of the canonical ensemble the 
physical temperature and pressure satisfying the zeroth law of thermodynamics were introduced. However, as it 
was shown in jMI^ HolllH such transformations of the variables including the entropy lead to the transition from 
the Tsallis thermostatistics to the extensive statistical mechanics of Gibbs or Reny i one. So such consideration 
does not pertain to the perception of the zeroth law of thermodynamics. In |l2| . a stronger assumption was 
evolved to make use of the nonextensive microscopic Hamiltonian of special kind like the nonextensive entropy 
form. In this case, the Hamiltonian depends on the temperature of the system that entails the changes of 
the thermodynamic relations and leads to the loss of self-consistency of the statistical mechanics. Closely 
connected with the problem of the zeroth law of thermodynamics is another one: the principle of additivity 
in the thermodynamic limit. In Abe [l^ . it was attempted to define the thermodynamic limit on a particular 
example of the perfect gas in the canonical ensemble. However, this limit was carried out incorrectly. The 
correct definition of the thermodynamic limit for the Tsallis thermostatistics in a particular case was given 



2 



in Botet et al. |lOl Il4| and for the general case was developed in [l5j. Note that an important criterion of 
correctness of the thermodynamic limit is the equivalence of all ensembles. Derivation of the fundamental 
equation of thermodynamics for the Tsallis thermostatistics on the base of the canonical distribution function 
was performed in p, TTl . In , it was proved that the microcanonical ensemble of the Tsallis statistical 
mechanics satisfies all requirements of the equilibrium thermodynamics if the entropic index l/{q ~ 1) is the 
extensive variable of state of the system. In the present Letter, we will show that similar results are carried out 
also for the canonical ensemble. 

The Letter is organized as follows. In the second section, the canonical ensemble and the derivation of the 
thermodynamic relations are given. In the third section, the performance of the thermodynamic principles in 
the thermodynamic limit on the example of the perfect gas is proved. 



II. CANONICAL ENSEMBLE 



Let us consider the equilibrium statistical ensemble of the classical dynamical systems of N particles at the 
constant temperature T, the volume V, and the thermodynamic coordinate 2: in a thermal contact with a heat 
bath. The system interacts weakly with its surroundings and only the energy can be transferred in and out of 
it. In order to determine the equilibrium distribution function, we consider the Tsallis equilibrium statistical 
entropy which is a function of the parameter q and a functional of the probing phase distribution function 
g{x,p): 

Q-Ql 



s = -k '^r, (1) 

. I - Q 



where dT = dxdp is an infinitesimal element of phase space, k is the Boltzmann constant and g e R is the real 
parameter taking values < 5 < 00. The phase distribution function is normalized to unity: 



gdV = 1. (2) 
In the classical statistical mechanics the expectation value of the Hamiltonian can be written as 

{H) = I qH dr. (3) 



The phase distribution function depends on the first additive constants of motion of the system. Nevertheless, 
the mechanical laws are not sufficient to determine it unambiguously. For this reason additional postulates of 
the equilibrium thermodynamics are required. To express the equilibrium phase distribution function from the 
macroscopic variables of state, we consider the thermodynamic method explored in ,1 5j . In the state of thermal 
equilibrium the macroscopic system is characterized by the fundamental equation of thermodynamics 

TdSth = dE+ pdV + Xdz - fidN, (4) 

where 5'th(r, V, z, N) is the thermodynamic entropy, z and V are the "thermodynamic coordinates" ; X and p 
are the associated "forces"; fi is the chemical potential and E is the thermodynamic energy of the system. In 
the canonical ensemble the fundamental equation of thermodynamics at the fixed values of T,V,z,N can be 
rewritten as 

(TdSth - dE)Ty...N = 0. (5) 

Then, to express the phase distribution function g{x,p) through the variables of state (T, V, z, N), let us replace 
in Eq. ij^l the equilibrium thermodynamic entropy Sth and energy E of the macroscopic system with the 
statistical ones ^ and J^Jl. Then, one finds 



T^dq+ I dV 



_ 5{H) 



5g 5g 



dg~^-^dH} = 0, (6) 



where the symbol d before the functions H, g and q is the total differential in variables (T,V, z, N). The 
statistical parameters q and g must be expressed through the variables of state of the system to provide the 
unambiguous conformity between statistical and thermodynamic entropies. Since dH = 0, dq — and dg = 0, 
we obtain 



(7) 
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where a is a certain constant. For the microcanonical ensemble it was stated in ref. jl5j that 

1 
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(8) 



where the parameter z takes the values — oo < z < — 1 for < g < 1 and < z < oo for 1 < g < oo and in 
the limiting case for q — \^ we have z — ±cx). Substituting Eq. J^l into Q and using again equation J^l to 
eliminate the parameter a we arrive at the following expression for the equilibrium phase distribution function: 



1 



1 A-if 



z + 1 kT 



where k = E — TS and it is determined from the normalization condition Q 

1 A- H] 



1 



z + 1 kT 



dT = 1. 



(9) 



(10) 



Thus, A is a function of the variables of state, A = A{T,V,z,N). For the Jaynes principle derivation of the 
phase distribution function (|5J|, see Appendix A. The expectation value (A) of the dynamical variable A{x,p) 
can be defined as follows (cf. ijSJl): 



A 



1 



1 A- H 



z + 1 kT 



dT. 



Using Eqs. © and l^, we can write the entropy (Q) as 

z (H)-A 



z + 1 T 

In the canonical ensemble it is convenient to introduce the free energy as a thermodynamic potential: 

(H) + zA 



F = E-TS 



z + 1 



(11) 



(12) 
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This is the Lcgendre transform of the energy with respect to the entropy of the system with dE/dS — T. 

At this point we have obtained a rigorous derivation of classical statistical mechanics. The quantum statistical 
mechanics is constructed in analogy with the classical one. In this respect, all classical notions are replaced 
by the quantum mechanical ones. A dynamic state of the quantum system is defined by a vector of state 
\^{t)), which is an element of the abstract Hilbert space £h- The dynamic variables are represented by linear 
hermitian operators A acting on the elements of the Hilbert space. In particular, the Hamiltonian H is the 
linear hermitian operator acting on the vectors of state |\E'(t)). In the quantum statistical mechanics the mixed 
states are considered. A macrostate thus appears as a set of possible microstates, which are set up by state 
vectors \^r{t)), r = 1,2, . . ., each with its own probability Wr for its occurrence, which are eigenvalues of the 
statistical operator, Q{t)\^r{t)) = Wr\9r{t))- The integral over phase space of the classical functions is replaced 
by the trace of the corresponding quantum operators. 

Let us show the connection between the statistical mechanics and equilibrium thermodynamics. In this 
respect, we derive the thermodynamic relations for the canonical ensemble from the general point of view. Let 
us consider now the classical case. Applying the total differential operator with respect to the ensemble variables 
(T, V, z, N) on the entropy ^ and the norm equation and using Eq. © one finds 



TdS = Xdz + / dgHdT 



where 



X^kT J g[l~ g^/^l-\ng^^')]dT. 



(14) 



(15) 



The repeated application of this differential operator to Eq. lO and substitution of the results into Eq. (|14|l 
leads to the formula 



TdS = d{H) - J gdHdT + Xdz. 



(16) 



By virtue of the parametrical dependence of the Hamilton function H on the variables V and N we have the 
fundamental equation of thermodynamics |^, lll| 



TdS = d{H) + pdV + Xdz - fidN, 



(17) 
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where 

Here, the property of the Hamilton function {dH/dT)v,z,N — {dH/dz)T.v,N = is used. So we have proved 
that the statistical entropy with the phase distribution function (|5J| completely satisfies the fundamental 
equation of thermodynamics Q which was a starting point of our derivations concerning the distribution 
function. In fact, the differential form (|17|l is sufhcient to prove the connection of the thermodynamics and 
statistical mechanics. 

The thermodynamic potential of the canonical ensemble (T, V, z, N) is the Helmholtz free energy Hence, 
the differential of F which follows from Eq. I|17|) can be written as 

dF = -SdT - pdV - Xdz + ^idN, (20) 

and the thermodynamic relations are 

^--m . .-m) . (-) 



The free energy as a thermodynamic potential can easily be calculated in the framework of the canonical 
ensemble and it is used in order to obtain the functions of the state. 

The fundamental equation of thermodynamics (|17|l provides the first and second principles of thermodynamics 

SQ = TdS, SQ = d{H) + pdV + Xdz - ^dN, (23) 

where 5Q is a heat transfer by the system to the environment during a quasistatic transition of the system from 
one equilibrium state to a nearby one. 

Let us now find the important quantity, the heat capacity, which is defined from the general rule, 5Q ~ CdT. 
According to the first and second laws of thermodynamics, in the canonical ensemble the heat capacity at the 
fixed values of T,V, z, N can be written as 

c.M.,v,.»,.(g)_^.(i^_ 



Note that all thermodynamic relations of classical statistical mechanics that we have just given are also preserved 
in the case of formalism starting from quantum mechanics. 

It is important to note that the first formalism of the Tsallis statistics 1] was rejected in the literature. 
For details, see Q. In the present study, we reconsider this initial formalism and show that the Tsallis phase 
distribution function was not accurately determined. It is claimed in 0] that the dependent function of the state 
(3(T,V, z, N), or the Lagrange multiplier, is an independent parameter of the distribution function. However, 
we show that the consistent derivation of the phase distribution function in terms of the variables of state 
{T,V, z, N), based on the Jaynes principle with the Lagrange function introduced in 0, leads exactly to our 
Eqs.Q and H1U|) . For details, see Appendix A. Note that in the limit, z ±oo, all expressions given above 
resulted in ones of the conventional Gibbs thermostatistics. 



III. THE THERMODYNAMIC LIMIT. THE PERFECT GAS 



The statistical mechanics must satisfy the requirements of the equilibrium thermodynamics in the thermody- 
namic limit when the number of particles of the system considered is very large, and the relative magnitude of 
surface effects becomes negligible. For thermodynamic systems, which are homogeneous at a macroscopic scale, 
the principle of additivity is valid: the various quantities of interest can be classified into either extensive or 
intensive ones under a division of the system into macroscopic parts. The extensive quantity, for instance, the 
entropy considered as a function of the extensive variables V, z and N is homogeneous of degree 1: 



S{T, XV, Xz, XN) = XS{T, V, z, N), 



(25) 
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where A is a certain constant. However, the intensive quantities, like the pressure p, the chemical potential fj,, 
and X, are the homogeneous functions of degree zero: 

fi(T,XV,Xz,XN) = fi{T,V,z,N). (26) 

Note that the temperature T is an intensive variable of state remaining invariant under such a subdivision 
of the system. This invariance property of the temperature guarantees the fulfillment of the zeroth law of 
thermodynamics. Here the thermodynamic limit denotes the limiting statistical procedure N oo, v = V/N = 
const, z = z/N = const with keeping the main asymptotic on TV. It is meant to make an expansion of the 
functions of the state in powers of the small parameter 1/iV {N ^ 1) with large finite values of the variables 
y, z. Then the extensive variables A can be written (a > 0) as 

AiT, V, z, N) |Ar^oo,„.i=const - N[aiT, V, ~z) + 0(Ar-")] = NaiT, v, ~z), (27) 

whereas the intensive variables 4> take the following form: 

(/.(T, V, z, N) U-.oo,.,z=const = HT, V, S) + OiN-") = 4>{T, V, z), (28) 

where v is the specific volume, z is the specific z, and a — A/N is the specific A. It is important to note that the 
limiting statistical procedures {V oo, z/V = const, N/V = const) and (z — > ±oo, V/z = const, N/z = const) 
are equivalent with the limit, N — > oo, given above. Note that in Abe 13J, the thermodynamic limit for the 
Tsallis statistics is not correct because the limits iV — > oo and |z| ^ oo are not coordinated among themselves. 
Note that after applying the thermodynamic limit to the functions of state the Boltzmann-Gibbs limit, z ±oo, 
is provided by expansion of these functions in powers of the small parameter 1/z holding only the zero term of 
the power expansion. 

In the canonical ensemble, how to prove from the general point of view the principle of additivity (see Eqs. (|25|l 
and (j2ni) and the zeroth law is not obvious. Therefore, we will illustrate explicitly the implementation of these 
principles on the foregoing example of the nonrelativistic ideal gas. In the framework of the nonrelativistic ideal 
gas of N identical particles in the canonical ensemble the functions of state can be explicitly expressed from the 
variables of state. Hence, the thermodynamic properties of the Tsallis statistics can thoroughly be investigated. 
In order to evaluate the expectation values of the dynamical variables, we use the method based on the integral 
representation of the Euler gamma function llll [la | . Let us investigate the thermodynamic properties of the 
nonrelativistic perfect gas of N identical particles in the thermodynamic limit (N —foo,v = const, z = const). 
The exact relations for it can be found in Appendix B. In the thermodynamic limit one easily confirms that the 
canonical partition function of the ideal gas for Gibbs statistics is simply 

Z'J''iT,V,N) (gve) (^)"' ^ Zg{T,v), (29) 

where Zq is the one-particle partition function, m is the particle mass and g is the spin degeneracy factor. Then 
in the thermodynamic limit and in the limit of Boltzmann-Gibbs statistics, z. — > ±oo, Eqs. (jBip and ljB2|) can 
be written as 



B{T,V,z,N) = (Zce^/^) 



B{T,vrz), (30) 



5|5^±oo = l--AnZG + 0{{l/zf), (31) 



where K/N — kTz{B~ 1) and z < —3/2 or z > with respect to the conditions of the integration method used. 
Note that the function of state B is intensive and A is extensive. The energy (IB3|) takes the following form: 



E{T,V,z,N) 3,^ B 
N ^ 2 

where e is the specific energy depending only on intensive variables and it is intensive. Consequently, the energy 
E is extensive (cf. Eas. l|27(l and (j^HJ)- Note that in the Boltzmann-Gibbs limit the energy (|32|l is reduced to 
E\i->±oa = {3/2)kTN. The free energy ||B4)| is found to be 



Eil2Xl^^kT~z{B-l) = f{T,vrz), (33) 

where / is the specific free energy which is intensive and F = A. In the Boltzmann-Gibbs limit the free 
energy H33fl takes the usual form F\~^^^ = —kThiZc. Then in the thermodynamic limit and in the limit of 



Boltzmann-Gibbs statistics the entropy IjBSp can be written as 



S{T,V,z,N) 
N 

•^Iz— ±00 



kz 
3 



B 



1 + ^ 



kN + k\n Zr 



Sc. 



(34) 
(35) 



where s is the specific entropy and Sg is the entropy of Gibbs statistics. The function s, however, is an intensive 
function of (T, w,z). Thus, the Tsalhs entropy (|34|) in the thermodynamic Umit is extensive. It is important to 
note that the Tsalhs entropy (|B5p for the finite values of iV and z does not satisfy the homogeneous condition 
(|25|l . for instance, in the case of A = 1/-^, and it is nonextensive. However, the Gibbs entropy, Sq-, is also 
nonextensive for the nonrelativistic perfect gas of iV identical particles in the canonical ensemble, because the 
canonical partition function Zq ^ (iV!)~^. See Appendix B. In the thermodynamic limit the Gibbs entropy 
(|^ of the perfect gas in contrast with its exact value is extensive, because Eq. (1^ has been implemented. 
Note that the values of the temperature and the specific volume are restricted, T > Tq and v > vq, hy the 
physical conditions S* > and Sq > 0, where the parameters (Tq, vq) are determined from the following equation 
ZG{To,vo)e^/^ = 1. 

Let us deduce the pressure p, the chemical potential fi, and the variable X in the thermodynamic limit. The 
pressure (|B6|I takes the form 



p{T,V,z,N) 



kT B 



" 1 + ^ 



= p(T,v,z), 



(36) 



where p is the intensive function and in the Gibbs limit we have pI^^^qq = kT /v. The function of state ljB7l 
can be written as 



X{T,V,z,N) = kT 



B 



J3_ 

2i 



1 - In ■ 



B 



_3_ 

22 



X(T,i;,5) 



and in the Gibbs limit it is reduced to X\^^j^^ 



0. The chemical potential l)B8|l is now 



m(T, V, z, iV) = kT- 



B 



_3_ 

22 



; In ■ 



B 



_3_ 

2z 



= KT,v,z), 



(37) 



(38) 



where fi is an intensive one and in the Gibbs limit we find that 

(jnU, and H36|l - H38|l in the thermodynamic limit yield the Euler theorem 

TS = E + pV + Xz- fiN. 



kT{l - InZc). Then, Eqs. 



(39) 



This equation enables us to interpret X as the "force" observable of the system conjugate to the "position" 
variable z. Moreover, Eq. H39|) allows us to write 



F = E-TS = -pV -Xz + ijlN. 
Its differential leads to the Gibbs-Duhem relation: 



(40) 



SdT = Vdp + zdX - Ndfi, (41) 

which shows that the variables T,p,X and fj, are not independent. The heat capacity of the perfect gas of N 
identical particles is deduced from Eq. I|24|l . Using Eq. H32(l we obtain 

C,,{T,v,S) = lkN (42) 

In the Boltzman-Gibbs limit the heat capacity takes its usual form Cvz\z^±oo ~ (3/2)fciV. 

To better understand the thermodynamic properties of the perfect gas, it is necessary to study the equilib- 
rium distribution of the momenta pi of the gas particles. The A'^-particle distribution function is defined in 
Appendix B. In the thermodynamic limit the single-particle distribution function (|B11|) takes the form 



7 



where rUeff is the effective particle mass written as 



B 

^ ' z 2 



gv 



kTe^ 



3/2' 



-1/z 



(44) 



In the Boltzmann-Gibbs Hmit, ^ef fl^^j.^ = ^n, and the single-particle distribution function is reduced to 
the Maxwell-Boltzmann distribution fip)\^^j.^ = (27rmA:T)~'^/^ exp(— p^/2mA:T). It is clearly seen that the 
effective mass (|44|l is rUeff > m for z < —3/2 and me// < m for z > 0. At z ^ —3/2 we have rUeff — > oo. Let 
us investigate the single-particle averages. The particle mean kinetic energy with the single-particle distribution 
function H43|) can be written as 

y \ ^ kT—^ = e. (45) 



2m / 2 m 

So the average kinetic energy is equivalent with the specific energy per particle H32|) . The average momentum 
of the particle and the highest probability momentum of the distribution f{p) can be written as 



(P) = J , = V2m,//fcT, (46) 



where p = \p\. Their ratio is {p)/php — '^/V^ ~ 1.13 as in the Gibbs statistics for which the average momentum 
and the highest probability momentum are {p)c = (8mfcT/7r)^/^ and {php)G — (2mA;T)^/^. The variance and 
the relative statistical fluctuations for the distribution (|43|) are 



{{Apr) = {p') {pY = kTm^ff 3 - - , (47) 



y/gA^ /3^ 



(48) 



or (5p K, 0.424. So from the investigation of the distribution function (|43() and its averages we arrive at the 
conclusion that in the Tsallis statistics the mean kinetic energy of the particles and the momentum are larger 
than their values in the Gibbs statistics, e > ea and {p) > {p)g, for z < —3/2 and smaller, e < ec and 
(p) ^ (p)gj for z > 0. Now we can give a physical interpretation for the variable of state z in the framework of 
the Tsallis thermostatistics if we consider the system of noninteracting particles as the system of the interacting 
quasiparticles with the effective mass We//- Then the total energy of the quasiparticle, e, is equal to the sum 
of the mean kinetic energy and the effective interaction energy Ae which can be written as 



/ p^ \ 3 , „ . 3 , „ rrie f f 



\ / = :;kT, Ae = -kT 

\2me/// 2 2 



1 



(49) 



So the effective interaction energy is positive Ae > 0, and the forces are repulsive for z < —3/2 and they are 
attractive Ae < for i > 0. Note that from Eq. I|49|l follows the physical interpretation for the temperature of 
the system T as the average kinetic energy of the quasiparticles with the effective mass m^f / . 

Fig. 1 shows the dependence of the single-particle distribution function (left panel) on the momentum p for 
the perfect gas of free quarks in the nonrelativistic approximation. The calculations are carried out for the 
system of quarks with two flavor and three color degrees of freedom, and the constituent mass m = 300 MeV 
at temperature T — 100 MeV and the specific volume v = 0.25/po, where po = 0.168/m~^. Note that the 
quarks at short distances are treated as almost free pointlike noninteracting particles because of the property 
of asymptotic freedom. It is remarkable that the single-particle distribution function at small values of the 
variable \z\ considerably differs from the limiting Maxwell-Boltzmann distribution. Such a behavior has really 
been caused by the sharp changes of the effective mass of quasiparticles in the dependence on z. This dependence 
can be seen even better in the right panel of Fig. 1 which shows the ratio of the effective mass to the constituent 
one vs. the variable z for two values of the temperature T. Figure 1 clearly shows that the variable z is the 
order parameter and the system is physically unstable in the region —3/2 < z < 0. For the microcanonical 
ensemble similar results were obtained in 15]. 

Our exploration of the perfect gas has brought the following points to the foreground: for the Tsallis statistics 
in the thermodynamic limit the principle of additivity and the zeroth law of thermodynamics are valid. The 
entropy (|34|l is a homogeneous function of the first order, it is an extensive variable satisfying the relation (|25|l 
with A = 1/A'' provided that the temperature T must be the intensive variable of state. It should, however, 
be noted that in the thermodynamic limit the equivalence of the canonical and microcanonical ensembles is 



io' 10' 10' 



p (MeV/c) 



-40 -20 



20 40 
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FIG. 1: The dependence of the single-particle distribution function on the momentum p (left) for the classical perfect 
gas of quarks with two flavor and three color degrees of freedom in the non-relativistic approximation for the different 
values of i = —2, —3, 3 and 2 (the curves 1, 2, 3 and 4, respectively) at the temperature T — 100 MeV and the specific 
volume V — 0.25/po- The ratio of the effective quark mass to the constituent one as a function of the specific z (right) for 
the values of T = 100 MeV (dashed line) and T = 200 MeV (solid one). The dotted lines correspond to the conventional 
Boltzmann-Gibbs thermostatistics. 



implemented. Giving constants of motion exactly or only as averages leads to the same results. For instance, if 
we express the temperature T through the variables {E, V, z, N) from Eq. as 



and substitute it into Eqs. H34|l and (|36() - (|38|l . then we obtain the results of the microcanonical ensemble derived 
in ref. |l5j| . It is important to note that in terms of the z variable and our thermodynamic limit the results for 
the perfect gas in the canonical ensemble of Abe et al. 0, ^| are the same as here. 



In this Letter, the canonical ensemble of the nonextensive thermostatistics introduced by C. Tsallis has been 
reconsidered. It is shown that the equilibrium statistical mechanics based on the nonadditive statistical entropy 
completely satisfies all requirements of the equilibrium thermodynamics in the thermodynamic limit. The unique 
non-Gibbs phase distribution function corresponding to the Tsallis entropy is obtained from the constraints 
imposed by the equilibrium thermodynamics laws. The microscopic foundation of the equilibrium statistical 
mechanics proceeds on the Gibbs idea of the statistical ensembles for the quantum and classical mechanics. The 
phase distribution function and the statistical operator depend only on the first additive constants of motion of 
the system. Here they were derived within a formalism based on the fundamental equation of thermodynamics 
and statistical definition of the functions of state. It allows us to avoid introduction of the controversial Lagrange 
multiplies. Nevertheless, it is shown that the distribution function derived from the Jaynes principle exactly 
coincides with ours if the Lagrange parameters are expressed through a set of independent variables of state of 
the system. The unambiguous connection of the statistical mechanics with thermodynamics is established. The 
equilibrium distribution function satisfies the fundamental equation of thermodynamics, the first and the second 
principles preserving the Legendre transformation. All thermodynamic relations relative to the Helmholtz free 
energy for the thermodynamic system in the thermostat are carried out. The heat capacity of the system 
was derived from the first and the second laws of thermodynamics. Note that in the fundamental equation of 
thermodynamics the new term related to the work of the conjugate force X at changes of the variable of state 
z appeared. In the limit z — s- ±00 the conventional Gibbs statistics is recovered. 

It is well known that the statistical mechanics should satisfy all requirements of the equilibrium thermody- 
namics only in the thermodynamic limit. Based on a particular example of the ideal gas we obviously proved 
the fulfillment of the principle of additivity and the zero law of thermodynamics for the Tsallis statistics in the 
thermodynamic limit. It was shown that all functions of state of the system are the homogeneous functions 




(50) 
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of the first degree, extensive, or the homogeneous functions of the zero degree, intensive. In particular, the 
temperature is an intensive variable and thus provides implementation of the zero law of thermodynamics. It 
should be marked that for the finite values of the number of particles N and the parameter z of the system 
both the Tsallis entropy and the Gibbs entropy are nonextensive functions of state, while in the thermodynamic 
limit they become extensive variables. The homogeneous properties of the functions of state allow us to find 
the Euler theorem and the Gibbs-Duhem relation. After applying the thermodynamic limit the expressions of 
the Gibbs statistics are obtained by the limiting procedure z — > ±00. The one-particle distribution function 
in the thermodynamic limit leads to the Maxwell-Boltzmann distribution function with the effective mass of 
particles me//. This distribution allows us to find the physical interpretation for the variable of state z as the 
order parameter of the interacting system of quasiparticles with mass mg/ / and the physical interpretation for 
the temperature T as the average kinetic energy of quasiparticles. The numerical example for the ideal gas 
of quarks with two flavor and three color degrees of freedom in the nonrelativistic approximation shows that 
the dense system of quarks in the dependence of values of the parameter z can pass from the strong coupled 
state to the repulsive state of quarks. In the framework of the ideal gas of identical particles the equivalence of 
the canonical and microcanonical ensembles in the thermodynamic limit is proved. This property is the key to 
identifying the self-consistency of the statistical mechanics and thermodynamics. 
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APPENDIX A: PHASE DISTRIBUTION FUNCTION 



Here, to derive the distribution function in the canonical ensemble we use the Jaynes principle 19]. In this 
respect, the Lagrange function can be written as 



Hq'] = ^ ~ " (/ q' dT-l^- [3 (^j g'H dT - {H)^ , 



(Al) 



where g' is the probing distribution function. After maximizing the Lagrange function (|A1(I . S^\^,^^ ~ 0, and 
using Eq. |^ to eliminate the parameter a, we arrive at the following expression for the equilibrium phase 
distribution function: 



(A2) 



where A = (H) —qS/k(3. Differentiating the function A and Eq. Q with respect to f3, and using the distribution 
function ljA2|l . one finds 



op op 



The parameter (3 can be related to the temperature 

1 _ dS dS/d(3 



T dE d{H)/dl3 
Then the distribution function (|A2|I takes the form 



kf3, 



l + {q-l) 



A-H 
kTq 



where A is determined from Eq. (j^J 



l + (g-l) 



A-H 



kTq 



kT' 



1 A-H 
z + 1 kT 



dT = 1. 



(A3) 



(A4) 



(A5) 



(A6) 



Note that Eqs. ljA5p and (jA6|l are identical with © and (|10|l . So the form of the distribution function in terms 
of the variables of state is independent of the method of derivation. 

Let us show that the distribution function expressed through the variables of state (T, V, z, N) in is 
equivalent to Eq. IjASp . So in the Lagrange function was written as 



a J g' dT- aP{q - 1) J g'H dT. 



(A7) 
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After maximizing (|A7|I we obtain 



g - - [1 - f3iq - 1)H] — , 
Z = I [I- l3{q^l)H]^ dT. 



(A8) 
(A9) 

To exp ress the Lagrange parameter f3 through the variables of state T,V,z,N, we use the method described 

^9-1 



) exp r 

in [l^l ■ FinaUy, for the Lagrange parameter f3 we get 



kTq 



(AlO) 



Substituting Eq. IjAlOj) into Eas. ljA8| l and (|A9|) and introducing the new function A in the following form: 



(All) 

we obtain the phase distribution function (|A5|) with the normalization condition (|A6|) . So the form of the 
Lagrange function does not disturb the distribution function in terms of the variables of state. 



APPENDIX B: THE FINITE PERFECT GAS 

Following the arguments given in Ref. we easily derive the norm function A from Eq. H1(J|) in the case of 
z < -1: 



1 + 



1 A 

z + lkf 



BiT,V,z,N), 



z < -1, 



(Bl) 



where Zc = ((qV)^ /N\)(mkT/2nh'^)^^^^ is the partition function of the conventional ideal gas of the 
Boltzmann-Gibbs statistics 0,123 and —z — > 0. In the case of z > 0, we obtain 



1 A 

z + lkf 



(z + i)^^r(z + i) 

T{z + l + lN) 



BiT,V,z,N), 



z > 0, 



(B2) 



where the new function B is introduced for convenience. In order to determine the energy of system, we insert 
the Hamilton function A(x,p) = H{x,p) into Ea. Hll(l and after performing integration, we obtain 



{H) = -kTN - ^ 
^ ' 2 1 



1 3 



N' 



(B3) 



where the function B is determined from Eqs. (|Bip and (|B2|1 . The norm function A and the energy {H) allow us 
to calculate the thermodynamic potential of the canonical ensemble, the free energy F. Substituting Eas. ljBl|) . 
(jB^ and (|nSl into Eq.|(^, we find 



F 



-kTz 



1 + -^^N 

' z+l 2^* 



Then, the entropy of the system can be easily obtained from Eqs.(J2J or I|1H|I : 

B 



S — kz 



1 



1 



1 3 
z + l 2 



N 



(B4) 



(B5) 



Let us calculate the pressure p, the chemical potential fi, and the variable X. Taking into account Eqs. (|21|l 
and IIB4II. one finds 



B 



_2 E 

z+l 2^^ ^ 



Differentiating Eq. (|B4|) with respect to z in conformity with Eq. (|22|l . we obtain 



X ^kT\l 



B 



1 + — -Af 

' z+l 2^* 



^z+l> 2^^ 

1 + — -A 

^ z+l 2^^ 



InB + Tp [a + j-N ] - ^(a) 



(B6) 



(B7) 
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where ipiv) is the psi-function which depends on arguments a = — z,7 = — 1 for z < — 1 and a = z+ l,7 = l 
for z > 0. Taking the derivative of Hli4|l with respect to N in conformity with Eq. (|22|l . we get 



/i — -kT 



B 



1 

z+l 



1 3^ -ln(7(z + l)B) + ^(^a + 7-7V 



2 + 1 2 



kT- 



B 



1 + 



z+l 2 



hi 



V 



mkT 



3/2' 



ip{N+l) 



(B8) 



The iV-particle distribution function of the classical ideal gas in the canonical ensemble for the Tsallis statistics 
can be written as 



/(pi, . . . ,pn) 



l^i=l 2m 

Z + l kT 



(B9) 



which is normalized to unity 



d^pi ■ --(fpNfipi, ■ ■ ■ ,Pn) = 1- 



(BIO) 



Then the reduced single-particle distribution function can be easily obtained by directly performing the integral 
over the momenta {p2, ■ ■ ■ ,Pn)- 



fiP) 
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(z + l)27rmkTB 



3/2 



1 



^ -|Z+|(Ar-l) 



1 2mkTB 



(Bll) 
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